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We study experimentally and theoretically structural defects which are formed during the tran- 
sition from a laser cooled cloud to a Coulomb crystal, consisting of tens of ions in a linear ra- 
diofrequency trap. We demonstrate the creation of predicted topological defects ('kinks') in purely 
two-dimensional crystals, and also find kinks which show novel dynamical features in a regime of 
parameters not considered before. The kinks are always observed at the centre of the trap, showing 
a large nonlinear localized excitation, and the probability of their occurrence surprisingly saturates 
at rsj 0.5. Simulations reveal a strong anharmonicity of the kink's internal mode of vibration, due to 
the kink's extension into three dimensions. As a consequence, the periodic Peierls-Nabarro potential 
experienced by a discrete kink becomes a globally confining potential, capable of trapping one cooled 
defect at the center of the crystal. 
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Stable, collective configurations that are non- 
perturbative, such as solitons, kinks, vortices, monopoles, 
and other structural defects, have been extensively stud- 
ied theoretically and experimentally in a variety of fields 
in physics, ranging from fluid mechanics, condensed mat- 
ter, atomic physics and optics, to high energy physics 
and cosmology [J-Q. These nonlinear solutions are lo- 
calized, particle-like objects that can propagate with- 
out dispersing, and acquire stability which can often be 
explained by their underlying topological nature. Such 
topological configurations can be formed during a first- 
or second-order phase transition 0-0 , for instance as in 
the case of vortex defects in a superfluid. Over the last 
decades, there has been significant interest in the study 
of topological defects in discrete systems, which can show 
similar nonlinear phenomena. The Frenkel-Kontorova 
(FK) model [7], for example, requires a periodic, one- 
dimensional lattice system, and features much of the ef- 
fects that are characteristic in its continuum counterpart 
(sine-Gordon field theory). In discrete systems, solitons 
are often referred to as kinks, translation-invariance is re- 
placed by lattice-periodicity and the kinks, propagating 
along the lattice, are subject to the so called Peierls- 
Nabarro (PN) potential [7[. The experimental study of 
discrete solitons has been demonstrated in several sys- 
tems, e.g., in waveguide arrays Q.and proposed with, 
e.g., cold atoms in optical lattices [9]. Since discrete de- 
fects possess substantially less degrees of freedom, they 
are in many cases experimentally favorable for precise 
control of their gap-separated, localized modes of oscilla- 
tion and in particular, for studying quantum coherence 
effects [loj and quantum information prospects lllll2|. 



In particular, cold trapped ions [13| offer a promis- 
ing platform for studying the formation and structure 
of topological defects, as well as for the exploration of 
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FIG. 1. Coulomb crystals, here with 31 ions (cj z /cj y ~ 1.34). 
(a) Top: a CCD image of flourescence light of individual ions 
trapped and laser cooled in a zigzag configuration. Bottom: 
circles indicate the fitted positions, while crosses depict the 
expected positions at the given trapping parameters, derived 
by MD simulations, the deviation remaining within the exper- 
imental resolution (see text), (b) As in (a), however, the ions 
in the center region forming an extended 'kink' defect. While 
the ions to the left of the center form a zigzag structure as in 
(a), the positions on the right are mirrored about the x-axis. 
Such a kink, linking two configurations of an odd number of 
ions in a circular trap, would be topologically protected. 



their quantum properties. When applying laser-cooling 
to trapped ions, they undergo a transition from a state 
of a chaotic cloud to an ordered structure, a Coulomb 
crystal [14| . Isolated from environmental disturbances in 
ultrahigh vacuum, such crystals forming one-dimensional 
chains permit investigating quantum effects at the fore- 
front of quantum metrology [151 ]. quantum computing 
[TR [13] and quantum simulations [l8l-[20| . By adiabat- 
ically altering the trapping parameters it is possible to 
generate phase transitions from the one- (lin ear) to two- 
zigzag) and three-dimensional crystals [2l|-[23j. Such 
structural transitions have been proposed to permit in- 
vestigating quantum phase transitions [23|, [24| . 
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The experimental creation of localized defects in ion 
crystals has been achieved during a crystallization of the 
ion cloud [l3[. It has been also suggested that a non- 
adiabatic change of the trapping potential could trigger 
the creation of defects during the phase-transition from 
one- to two-dimensional crystals 25, 26|. Realizations of 
the FK model with ion strings have been suggested in 
27-29|. The quantum mechanical properties of kinks in 
ion traps have been theoretically explored in [25| , and the 
connection with the inhomogenous Kibble- Zurek mecha- 
nism in 26|, 30|. Currently the latter is under investiga- 
tion in several groups, e. g. (3l|, having created kinks by 
quenching the radial confinement of the trap. 

In this Letter we present the experimental realization 
of predicted, extended 2D-defects as well as a new class of 
quasi-3D defects, both induced by laser cooling [32| and 
long lived in the Doppler-cooling regime. In the ther- 
modynamic limit, as well as in circular configurations 
(25I . [33|, similar discrete solitons (kinks) become topo- 
logically protected. With respect to these we regard the 
observed configurations as 'topological' defects. For the 
3D-kink we directly observe a highly nonlinear oscilla- 
tion at its localized, gapped mode and reveal that the 
probability to observe a kink at equilibrium shows an in- 
triguing saturation at ~ 0.5. We study numerically the 
structural and dynamical properties of these kinks and 
find that their extention into 3D causes a modification 
of the Peierls-Nabarro potential, featuring an inherent 
trapping of the kink at the centre of the Coulomb crys- 
tal, and resulting in the observed saturation in the kink's 
occurrence. 

The experimental setup consists of a linear Paul trap 
(fi rf 2tt • 6.22 MHz) [3, EH trapping Coulomb crys- 
tals of different size of choice (here 10 - 65 photoionized 
24 Mg + ions). Experiments were carried out at a single ion 
trapping frequency of uj x ~ 2tt • 56kHz (axial direction), 
while the radial frequencies were varied in the range of 
uj y « 27r-320kHz to 27r-630kHz, corresponding to a total 
depth of the rf-trapping potential energy of the order of 
ke 10 4 K, ke being the Boltzmannfactor. The radial in- 
degeneracy uj z /uj y is tuned between 1 and 1.4. The ions 
are Doppler cooled via one laser beam, slightly tilted (~ 
5°) with respect to x-axis. It drives the S1/2-P3/2 tran- 
sition (A= 280 nm, natural linewidth T « 2tt-42MHz) 
and is detuned byA^2-T/2ata saturation of about 
0.2 (Doppler cooling limit T D « 1 mK at A = T/2). 
Data is acquired with a CCD-camera detecting fluores- 
cence light of 24 Mg + to determine the ion positions (see 
figure [T]). The resolution is given by the magnified (10 x) 
pixelsize (0.8 /im) of the image on the CCD-camera. In- 
dividual images are integrated for 200 ms. By increasing 
the axial confinement (decreasing the radial) we repro- 
duce increasingly complex structures, from a linear chain 
to 3D crystals. At constant confinement, the structure 
and its dimensionality depend on the number of trapped 
ions, due to the impact of their space charge density on 
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FIG. 2. Coulomb Crystals containing 50 ions, similar to those 
depicted in figure [TJ however for reduced asymmetry of the ra- 
dial confinement, cj z /cj y w 1.05. (a) CCD images of two crys- 
tals obtained for identical experimental parameters. Top: a 
pure zigzag featuring the minimal energy configuration. Bot- 
tom: two central ions clearly deviate and show a blurred ex- 
tension in the radial direction. As in figure[T]3, the right parts 
of the crystals are close to identical while the left part of the 
lower crystal depicts a flipped 'zagzig' structure, separated by 
the two central ions, (b) Integrated numerical results for the 
structure depicted in (a)-Bottom, projected onto the crystal 
plane. The agreement of theory and experiment remains com- 
parable to figure [T] (c) Numerically derived components of 
ion oscillations in the localized, low- frequency normal mode 
of the kink (see figure [4)3) in (perpendicular to) the crystal 
plane are depicted on the left (right). 



the total potential. 

To calibrate our system and to gain a deeper under- 
standing of the structure and the dynamics within the 
crystals, we numerically simulate our experimental re- 
sults by a Molecular Dynamics (MD) code [36]. We take 
into account the full time-dependent trapping potential, 
and run an optimization routine which considers all trap- 
ping parameters and the projection of the crystal plane 
(xy) on the CCD-chip chosen roughly perpendicular to 
the z-axis. We reproduce all observed crystals consis- 
tently, with a maximal average deviation of 0.5/im per 
ion, which is less than the pixel size. We obtain identical 
results in a simulation using a constant harmonic poten- 
tial, verifying that micromotion [37} does not alter our 
findings. 

In the investigated regime of trapping parameters, for 
a number of ions between 20 and 60, zigzag structures 
with and without local defects occur (see figure [I]). The 
duration after which a zigzag-structure gets destroyed, 
for example, by a collision with constituents of the back- 
ground gas, amounts to 10 seconds. The lifetime of struc- 
tural defects within crystals amounts to the same order 
of magnitude and exceeds the natural timescale for os- 
cillations of the normal modes in the Coulomb crystals 
(~ 1/cj x ) by five orders of magnitude. The two differ- 
ent structures imaged by the CCD-camera, depicted in 
figure [H are obtained for uj z /uj y ~ 1.34, where uj z is the 
in-plane and uj y is the out-of-plane radial trapping fre- 
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FIG. 3. Probability for the occurence of a pure zigzag and 
a zigzag with a kink as depicted in figure [2] formed during 
a nonadiabatic, cooling-induced transition from a cloud to 
a crystal, for identical trapping parameters, plotted in de- 
pendence of the number of ions. With 27 ions and more the 
minimum energy configuration at our set of trapping parame- 
ters is a zigzag structure (blue circles). Starting with 46 ions, 
we observe an increasing probability for a zigzag structure 
with one kink (red rhomb) at the centre of the crystal. No 
simultaneous occurrence of two or more defects is observed. 
Surprisingly, the probability for observing one kink saturates 
at rsj 0.5 (dashed line) for up to 56 ions (see explanation in 
the text). Error bars represent the statistical error (lcr). 



quency. Whereas figure QJi presents the established zigzag 
configuration, which is the global minimum of potential 
energy, figure {\]p depicts a local minimum, incorporating 
an extended kink, as predicted in [25j. The zigzag struc- 
tures considered so far are extended in two dimensions 
(xy) only, independent of the presence of a kink. Keeping 
the axial confinement constant and lowering cj z /u; y , we 
create qualitatively new defects (figure^). For the rest 
of this paper we discuss their occurrence and properties 
on the exemplary set of parameters of oj z /uj y « 1.05. In 
this regime, zigzag structures without a kink occur start- 
ing with 27 ions and remain two-dimensional with up to 
52 ions. With 53 ions, a structural phase transition into 
3D can be deduced from the numerical results, the space 
charge of the ions defocussing a few central ions slightly 
out of the zigzag plane (by <l/im). Further increasing 
the number of ions causes a growing fraction of the ions to 
extend into the z-direction, however, the global minimum 
configuration (which remains close to a planar zigzag) is 
still unambigously identifiable. For more than 56 ions no 
zigzag structures are observed anymore, and we observe 
projections of complex, helical structures. 

Starting with 46 ions, we detect the random occur- 
rence of local defects mainly involving two ions. The di- 
rectly observable features of the novel class of kinks are: 
(i) They sustain a large motional amplitude of the two 
ions along the y-axis, comparable to the radial extension 
of the crystal, causing the blurring of their images (see 
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FIG. 4. Dispersion relations for a crystal of 50 ions as in fig- 
ure O with (red rhombs) and without (blue circles) a kink, 
(a) The spectra of normalized frequencies cj n /cj x of the 150 
normal modes remain close to identical for the most part, (b) 
A magnified inset showing that the kink configuration has a 
low- frequency mode separated by a gap. It is localized to the 
central ions only, and causes their blurring in the xy-plane 
(see figure [2ji). (c) The normalized frequency of the local- 
ized low-frequency mode cji ow for a given number of ions, can 
be controlled by tuning the ratio uj z /uj y (here reaching up to 
cjiow ~ 2cj x )- When the frequency of the lowest mode ap- 
proaches zero, the configuration including the kink undergoes 
a structural phase transition. 



figure 2a). (ii) The probability for the occurrence of a 
kink increases with increasing the number of ions, and 
saturates at ~ 0.5 for more than 50 ions (see figure [3}. 
(iii) Kinks occur at the center of the crystal only. In 
order to elucidate feature (i), we use MD simulations as 
described above and reproduce the positions of the ions 
and their dynamics in the xy-plane (see figure Wp)- The 
crystals including one defect remain nearly identical to 
a pure zigzag, except at the center. There, despite the 
stiffer uj z > uj y confinement, the central ions extend into 
the transverse (z) direction (see figure^), even for pa- 
rameters such that the zigzag crystal without a defect 
is still planar. We further investigate the properties of 
the kinks numerically applying two separate approaches: 
(1) a linearization of the oscillations of the ions at small 
amplitudes; (2) simulating the adiabatic dynamics of the 
kinks as they move along the lattice. 

The dispersion relations for the pure zigzag and the 
crystal with the centered kink (figure [2} are shown in 
figure [H For the largest fraction, both structures ex- 
hibit indistinguishable mode spectra. However, whereas 
the lowest mode of the zigzag is the axial center of mass 
mode (with the frequency cj x ), the kink features a gapped 
mode below cj x (see figure [4^, b), localized on the cen- 
ter ions. A dynamical simulation assuming a thermal 
energy of kpTo for all modes and integrating a 'fluores- 
cence' image from the numerically obtained ion trajec- 
tories, reproduces accurately the blurring of the central 
ions (see figure [2J3). We numerically confirm that this 
effect is not caused by micromotion or by the near-axial 
incidence of the cooling laser, and that it remains un- 
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affected assuming a realistic equilibrium temperature of 
~ 5Td- In fact, the apparent blurring of the center ions 
is reproduced exactly, even assuming an excitation of the 
low-frequency mode only. A careful analysis of the dy- 
namics [36| reveals that the large-amplitude motion is 
due to the highly anharmonic oscillation of the ions on 
the localized low-frequency mode of the kink. In addi- 
tion, the frequency of the low energy mode, cji ow , can 
be controlled by tuning the radial trapping frequencies 
u z /ujy over a range of to ~ 2u; x (see figure Hh). 

To further elucidate features (ii) and (iii) of the ob- 
served kinks, we focus on this mode that describes the 
ion oscillations near a minimum of the Peierls-Nabarro 
(PN) potential. For ion cr ysta ls of constant ion density, 
e. g . in a circular trap 25|, [33| or in an anharmonic trap 
[38|, the PN potential is periodic with the inter- ion dis- 
tance. For a zigzag in a harmonic trap, the space charge 
density and the radial displacement decrease out of the 
centre of the crystal. If uj z /uj y ^> 1, the depth of the 
local minima of the PN potential decreases progressively 
towards the sides [26|, [30|, where the kink can dissolve 
into the linear chain at each of the ends, after a random 

walk jasl. 

To analyze the motion of the 3D-defect along the lat- 
tice for our case of uj z /uj y > 1, we run a dynamical MD 
simulation for a kink created initially at an off-centered 
position. We follow its adiabatic dynamics while cooling 
the motional degrees of freedom of the ions to overdamp 
the dynamics. We find that the kink 'slides' down to- 
wards the center of the crystal within an effective trap- 
ping potential. Tracing the 'downhill' motion of the kink 
allows to calculate the potential along the adiabatic kink 
trajectory. We reveal that the effective PN potential is 
no longer periodic, but rather modified to an overall har- 
monic potential of a global depth quadratically depend- 
ing on the total number of ions (see figure [5b). Similarly, 
since the planar extended kink of figure [T] expands more 
at its centre than the zigzag configuration, the PN po- 
tential is modified in the crystal plane such that this kink 
is stable only at the centre. 

For 44 ions, we derive from the simulation the global 
potential depth for a kink within the crystal to be roughly 
lOkeTD (figure EK) . This is the height of the barrier for 
the kink's escape to the sides where it is lost, and the rel- 
evant degree of freedom for the propagation of the kink 
is its localized mode. The minimal equilibrium temper- 
ature of the ions in the experiment can be estimated to 
amount to a few times Td- Therefore, during crystalliza- 
tion and formation of the defects, a kink can be created 
with sufficient energy to escape its trapping potential. 
However, the potential depth for a kink rises with the 
number of ions such that for 46 ions, it already amounts 
to ~20kBTD, permitting to trap the defect with higher 
probability within the Coulomb crystal. Since the mini- 
mum configuration has a Z2 symmetry (of a zigzag and 
its mirror-image), the defect following a kink must be 
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FIG. 5. Simulation of the effective PN potential energy of 
a kink (see figure [2}, (a) dependent on its distance from the 
center of the Coulomb crystal. The potential is symmetric 
about the trap center, and reaches its maximum roughly at 
the final values shown (dependent on the total ion number). 
Beyond the related displacements, the kink will 'slide' to the 
side and dissolve, (b) The total potential barrier for the kink 
to escape its trap by reaching the sides is given in dependence 
of the number of ions, fitted by a quadratic function (blue 
line). 



its anti-kink. All defects created within the trapping re- 
gion will slide to the center and pairwise annihilate, ei- 
ther leaving one centered kink (if the initial number of 
kinks was odd), or none. Although we are not yet capa- 
ble to observe this process directly, these considerations 
provide a consistent explanation for the experimental oc- 
currence of one centred kink remaining at ~ 0.5 for the 
whole range of 50-56 ions (see figure[3]), with the small no- 
ticeable oscillation possibly indicative of the underlying 
dynamics. Additionally, it will be interesting to analyze 
the interplay of thermal effects with the increase in the 
number of kinks being formed during the crystallization 
as the number of ions is increased. 

Entering the quantum regime for large Coulomb crys- 
tals and 2D structures by exploiting the kinks' localized 
properties is another interesting step. Cooling the lo- 
calized gapped mode to the motional ground state and 
studying its coupling to the bath consisting of the rest 
of the modes would permit to explore mesoscopic open 
quantum systems. The controllable trapping parameters 
offer the possibility to shape the characteristics of the 
kinks, e.g. by tuning the frequency and gap of the lo- 
calized mode as described above. The 'kink trap' storing 
kinks at the center of the Coulomb crystal will be ex- 
ploited for studying the interaction of different classes of 
defects (such as a kink combined with a mass defect or 
a second kink protected by a mass-defect, that we have 
observed and will describe elsewhere [36| ) . These cannot 
annihilate, but yet are forced into a short range interac- 
tion in a quasi-deterministic way. 
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